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ABSTRACT 

The problem of determining the trajectory geometry which minimizes the propellant expendi- 
ture necessary to attain specified end conditions, for the planar ascent of a constant thrust 
rocket-powered vehicle in vacuo , is treated by the second order Euler equations of variational 
calculus. 

A series of equations is derived which allows a rapid isolation of the initial values of 
certain Lagrangian multipliers (which are necessary to attain prescribed end conditions) via an 
approximation using a parabolic turn program. 

A restricted (though not approximate) solution of the general result is derived which does 
not involve these multipliers. This solution is of a simple nature, and may find application to 
"on-board M adaptive guidance computers. 
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Integration Constants 

Fundamental Function Appearing in the Euler Equations 
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The Quantity to be Minimized (Maximized) Including all Applicable Constraints 
Constants Appearing in the Parabolic Turn Program 
Planetary Mass 
Vehicle M ass 

Magnitude of the Radius Vector 
Magnitude of the Thrust Vector 
Time 

Magnitude of the Velocity* Vector 

Any of the Constrained Variables Appearing in F 

The Control (State) Variable jS 

Angle of Attack; Measured from the Velocity Vector to the Thrust Vector 

The Sum of the Angle of Attack and the Flight Path Angle; the Thrust 
Orientation Angle Measured from the Vertical to the Thrust Vector 

The Flight Path Angle; Measured from the Vertical to the Velocity Vector 
Time Dependent Lagrange Multipliers 

Central Angle (Longitude) Measured from Launch to the Local Radius Vector 
Denotes a Final (Cutoff) Point 

Denotes an Initial Point. Usually First Stage (Booster) Cutoff. 

Denotes the Maximum Value of a Quantity 


Time derivatives are denoted by the dot notation. 
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SUMMARY 

The problem of determining the trajectory geometry which minimizes the propellant expendi- 
ture necessary to attain specified end conditions, for the planar ascent of a constant thrust 
rocket-powered vehicle in vacuo t is treated by the second order Euler-Lagrange equations of 
variational calculus. 

The normal procedure of eliminating the time dependent turn program as a function of certain 
Lagrangian multipliers is reversed in that the multipliers are expressed in terms of the turn pro- 
gram and its derivatives. This allows an isolation of the (generally unknow r n) initial values of 
the multipliers in terms of a functional approximation of the optimal turn program. The approxi- 
mation presented (which is certainly neither unique or sacrosanct) assumes that the angle of 
attack is a parabolic function of time. 

This procedure yields another valuable result. A second order differential equation govern- 
ing the thrust orientation angle as a function of time which contains no Lagrangian multipliers 
is derived. This solution, a special case of the first result, contains the thrust orientation angle 
and its first derivative as the only variables having unknown initial values. These quantities 
are generally quite easily approximated by an experienced worker. Few people have an intui- 
tive feel for the initial values of the Lagrange multipliers. 

These solutions are described from the "degrees-of-freedom’ ’ viewpoint and numerical 
applications are indicated. 

SECTION I. INTRODUCTION 

Since the problem of optimal thrust programming to obtain maximum altitude for a given 
rocket powered vehicle was first po'sed by Goddard, over three hundred studies have been devoted 
to optimal rocket trajectories. These studies have considered optimal thrust programming and 
optimal trajectory geometry to obtain maximum or minimum values of some "payoff” at the end 
point. 

This "payoff” may be almost any desired function of the end conditions such as maximum 
range (for ballistic missiles), maximum altitude (Goddard’s problem for sounding rockets), maxi- 
mum payload (for satellite carriers), maximum endurance (for rocket boosted gliders), etc. Most 
authors have restricted the treatment of such trajectories to two dimensions, usually commenting 
that an extension to three dimensions is a simple matter (this is true in theory, but less straight- 
forward in practice). The problem of including drag effects has been treated extensively and a 
number of analytical drag models have been investigated. Finally, almost all writers in this 
field have considered simultaneous optimization of the trajectory geometry and burning program. 

The situation is modified considerably if operating computational programs are considered 
rather than theoretical studies. These programs are usually of a preliminary design nature since 



the general problem has not been formulated to a degree which allows inclusion of numerous 
important effects. On this level, the restriction to planar flight is a reasonable assumption, 
especially if we consider only upper stage trajectories. Drag is usually not considered for sev- 
eral reasons. First of all, this would require inequality constraints on the angle of attack, dyna- 
mic pressure, angle of attack at cutoff, etc. Another restriction is that analytical drag models 
tend to be rather crude while table look-up is difficult to apply to variational formulations so that 
the actual procedure for introducing drag into the calculations is not always clear cut. Another 
point of departure between theory and practice is that numerical calculations seldom are involved 
with variable thrust along the trajectory. This specialization is due to the fact that most exist- 
ing engines have a constant thrust level at the present day state-of-the-art. 

Even under the above restrictions there are few calculus of variations trajectory shaping 
procedures actually used. This is principally because the operation of such a program requires 
that the initial values of certain Lagrange multipliers (which are used in the formulation of the 
problem) be known. The general procedure is to approximate these initial values and conduct 
numerical iterations until the correct initial values are obtained. This approximating technique 
is usually effective only if the worker is quite skilled from past experience. The length of time 
required to attain this ability is usually rather long. 

This paper deals with a method of determining the initial values of the Lagrange multipliers 
without previous experience by any person familiar with a reasonable trajectory computation pro- 
cedure.* Furthermore, another optimal turn program which does not involve these multipliers is 
derived as an important by-product of the formulation. This turn program is more general than 
the so-called bilinear tangent solution in that the assumption of a flat earth is not invoked in the 
present development. 

The author wishes to express thanks to S. Peter Gary without whom this paper would never 
have been started, and Larry G. Singleton without whom it could never have been completed. 


SECTION IL EQUATIONS OF MOTION 

Let r and \fj be the polar coordinates which describe the position of the center of mass of a 
rocket-powered vehicle above the (spherical) earth. Let V be the magnitude of the velocity 
vector and -9- the angle from a unit vector in the radial direction to a unit vector in the direction 
of the velocity (the so-called flight path angle). Similarly, define a (the angle of attack) as the 
angle from the velocity vector to the thrust vector. Denoting the magnitude of the thrust by T, 
the mass of the vehicle by m, and the gravitational parameter of the earth by GM, the equations of 
motion, for vacuum flight, are: 


• T GM 

V = — cos a 7 T 


( 1 ) 


& 


T 


mV ’ V 

r — V cos & 

; v . 

if/ = — sin & 
r 


sin & / GM V 2 

a + — — ( — “ — 




( 2 ) 

( 3 ) 

( 4 ) 


* The paper is generally orientated toward maximization of payload for a satellite carrier, but 
has other applications as well. 


2 



Most variational formulations proceed directly from these equations; however, it will be 
convenient to modify this approach for our purposes. 

Taking the ratio of equation (4) to equation (3) gives 


Squaring equations (3) and (4), and adding the results we find 

V =(? + r 2 i 2 )' A «*) 

Now differentiating equation (3) 

zJL = r = V cos S- - V » sin » (7) 

dt 


The V and $ terms can now be eliminated from equation (7) via equations (1) and (2) which results 
in 


T 


T . GM 

— cos (a + -9*)- — 2 + 
m r 



9 


( 8 ) 


As a final modification of the last equation, we now eliminate V and sin 9 by use of equations 
(5) and (6). Thus 


T v GM ; 2 

r = cos (a + 9 ) - — j- + t xp 

m t 


(9) 


A similar differentiation of equation (4), followed by elimination of V y V, 9, , and 9 , gives 

T , v 2 r xb 

if; = sin (a + 9 ) - 


( 10 ) 


It is interesting to note that the equations of motion in vehicle-fixed coordinates[ equations 
(1) through (4) ] convert to polar coordinates form [equations (9) and (10)] by simple differentia- 
tions. 

In order to simplify the following algebra, we now substitute 


(3 - 9 + a 


into equations (9) and (10) obtaining 


]_ 

m 


$ 


cos - 

GM 

— 2 r + r 
r 

sin ft - 

2 r xli 
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;.2 


(ii) 


( 12 ) 

( 13 ) 
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Our variational formulation will assume equations (12) and (13) as a starting point. 


SECTION III. VARIATIONAL FORMULATION 

The problem at hand is to determine (3 as a function of time such that the "payoff” (the 
quantity to be maximized) is a maximum (or, more correctly, stationary) value. The payload* 
will be defined to be our "payoff”. 

As is well known, maximizing the payload is equivalent to minimizing the propellant expend- 
iture required to attain orbit. Denoting a final point with an / subscript, an initial point with a 
0 subscript, and the expended propellant by m we have 


= m Q - m { = m Q - (m Q - m t { ) = m t f 


(14) 


where t is the time and m is the (constant) mass flow rate. 

Our problem is now seen to be of the Mayer type; that is, we minimize t f subject only to 
the constraints of equations (12) and (13). 

Defining 




o GM } 2l 
cos j3 + — - r tA 


+ A, 


ilf sin R + 

m r r r 


2r xjj 


= F [ r, r, r, ifj , \jj , t ] = 0 


(15) 


where Aj and A 2 are time-dependent Lagrange multipliers, we may now define the problem by 
seeking the stationary values of 

H = m t f + f F [r (t), r(t), r(t), ijj (t)> ijj(t), fi(t) 4 t] dt (16) 

J *0 

Application of the Mayer formulation to the above expression shows that among the neces- 
sary conditions are the following Euler equations (Ref. 1). 


d 2 ( d F 
dt 2 \d x 


d__ (d_F\ d_F_ _ 
dt \d x, J + d x . 


(17) 


where x { = (r, ip ) and 


d F 
d y j 


(18) 


where y is a state variable, which for our case is ^3. 

From equation (15), it may be immediately noted that \p does not appear, and we thus obtain 
a first integral of the form 


d 

Tt 



d F 
d ip 


+ C, 


= 0 


(19) 


* The term "payload” is used to denote gross weight in orbit. Thus our definition lumps to- 
gether the structure of the last stage, residual propellant, actual payload, etc., under the general 
heading of payload. 
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where C x is a constant of integration. 

Applying equation (18) to equation (15) with y i - ft we find 

— ( \ , sin fi - cos fi) = 0 (20) 

m \ 1 r ' 

T 

For * 0 and cos 6 * 0, equation (20) now gives 

m 


\ 2 = r tan £ 


( 21 ) 


The variational equation for the xfj coordinate may be derived from equations (15) and (19). 
For the first term 


d (d F 


dt \d if/ 


= A, 


( 22 ) 


and, for the second term 


JJ~ = " 2 A ‘ 


Inserting equations (22) and (23) into equation (19) we have 

A 2 +2 Aj f ^ ~2 A 2 ~ + C j = 0 

Now turning to equations (15) and (17) for x i = r, we find 


(23) 


(24) 


d 2 



(25) 


dt 






+ 2 



2 A 2 r j/r 
r 2 


( 26 ) 


<9 F 

(9 r 



A x i/f 2 + A 2 -^2 s ^ n - 2 A 2 




(27) 


Inserting equations (25), (26), and (27) into equation (17) gives 

Ai “-2A 2 ^3“2A2-^-“-2 Ai ^ - Ai \f/ 2 + A 2 ~2 s * n = 0 (28) 

Equation (28) may be rewritten by replacing A 2 by the expression obtained from equation 
(24), if/ by the right-hand side of equation (13) and A 2 by the right-hand side of equation (21). 
This results in 


Ai + Aj 



sin 


j8 tan jS + 2 Cj 




= 0 


(29) 
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Equations (12) and (13), the equations of motion, along with equations (24) and (29), govern- 
ing the Lagrange multipliers, and equation (21), which relates the steering program to the La- 
grange multipliers, are five equations in the five unknowns r, ijj, , A 2 , and These equa- 
tions, along with the initial values r 0 , r 0 , i// 0 , j}j 0 , /3 0 (which are known), (A^o , (Ai) 0 , (A 2 )o> 
(which are generally unknown) and C i (which is constant though unknown) are sufficient for a 
complete solution of the problem. 

SECTION IV. FURTHER DEVELOPMENT OF THE TURN PROGRAM EQUATION 

We shall now develop relationships between the Lagrange multipliers and the angle /3 
which permit an estimate of the initial values of the multipliers in terms of an approximate turn 
program . 

Taking the time derivative of equation (21) gives 

A 2 = A x t tan /3 + A x r tan /3 + A j r sec 2 /3 (30) 

Eliminating A 2 from equation (30) via equation (24) results in 

A i r tan jS-Aj r tan /3 + Ai r /3 sec 2 £J+2A 1 ri^+C 1 = 0 (31) 

Differentiating equation (31) now produces 

t tan fl + 2 A x r (0 + /3 sec 2 /3) 

+ A j [ — r tan /3+2ri/f+2r^+r(/3 sec- 2 + 2 ^ 2 sec2 /3 tan /3) ] = 0 (32) 

Aj , A j , V, and If may now be eliminated from equation (32) through the use of equations (29), 
(31), (12), and (13), respectively. Substitution from these equations and algebraic manipulation 
of the resulting form give 

cos2/3 | A x 3 sin 2 /3 - 4 r if/ 2 + — sin /3 tan /3 

+ r ( f3 tan (3 - 2 ft 2 - 6/3 if/ ) + 2r/3 tan 

- 2 Cj (t/y + /S)| =0 (33) 

The possible solutions and implications of equation (33) will now be discussed. 

SECTION V. TURN PROGRAMS 

Equations (12), (13), (21), (24), and (29) correspond to the results obtained when variational 
methods are applied to the study of rocket-powered vehicle trajectories. Equation (33) is a con- 
verse solution, insofar as the multipliers A j and C j (or at least their ratio) are determined, if 
the turn program is known as a function of time. As will now be shown, a combination of the two 
points of view yields more information than either approach produces separately. 

Several possible solutions to equation (33) immediately present themselves. The first of 
these, found by setting cos /3 = 0, was formally excluded by the derivation of equation (21). We 
proceed directly to the following possibility: 

[ CM 2 * 2 T • ? • • 

3 ~pr s * n P ~ 4 T \jj + — sin j3 tan /3 + r ( /3 tan /S-2/3 — 6 /3 v> ) 

+ 2 r /3 tan (3 — 2 Cj + /3) = 0 (34) 
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OPTION ONE 


The most general case for a solution of this equation assumes that 

Aj y 0 


( 35 ) 


and 


C i * 0 


(36) 


If equations (35) and (36) are fulfilled, we now have a method of determining all initial values of 
the Lagrange multipliers. 

To illustrate, the pertinent equations of Sections III and IVwill be gathered together. These 

are: 


r = 

T c 

— cos t 
m 

. CM 

r 2 

- + r i/r 

2 




4> = 

T 

mr 

sin 

P- 

2 r tj/ 
r 





A t + 

K 1 

(3 rjt 

2 -2 

CM \ 
r 3 ) 

-K 

— sin p 
mr 

tan jS + 2 C j 

0 

II 

A 2 + 

2 A, 

r* 

- 2 A 2 

T 

+ 

r 

Cl ■ 

= 0 



*2 = 

*i 

r tan 

P 






A, = 

A, 

r 

LT ‘ 

- /3 sec 

CSC 

p- 

2 \fj cot jS 

- cotfi 

r ' 



m = m 0 — m t 


( 12 ) 

(13) 

(29) 

(24) 

( 21 ) 


(37) 

(38) 


where equation (37) was derived from equation (31) and equation (38) gives an explicit relation- 
ship between the vehicle mass and time.* 

Once the initial values of the physical variables r, r, ijj, ijj, T and m (or t) are specified, 
the only available choice of initial conditions to obtain desired end conditions are (A^q, (Ai) 0 . 
(A 2 )„, C l > and the burnin g time h~ l o • The choice of these parameters then allows a numerical 
integration of equations (12), (13), (21), (24), (29), and (38) to be carried out. One of these values 
is arbitrary, however, as can be seen by noticing that equations (21), (24), and (29) are homo- 
geneous in the multipliers (including C, ). Thus, we are left with one arbitrary choice, three 
unknown choices, and the burning time (which may conveniently be used to specify energy or 
velocity). The choice of which initial value is arbitrary is actually more important than what 
value is chosen. Numerical work performed with this set of equations indicates that (Ai) 0 is a 
good choice since the end conditions are rather insensitive to variations in this parameter (with 
another initial value fixed).** 

* It is unnecessary to include equation (38) as a constraint in the formulation of equation (15) 
since we assume an explicit time dependence in the F function. 

** A good choice for upper stage trajectories of large modem vehicles is (A 1 ) 0 = 50.0. 
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Proceeding on the ^assumption that \ x is arbitrarily fixed, we now further assume that the 
initial values of fi Q , /3 0 , and /S 0 are known at least approximately. From equation (34) we may 
now determine the value of C \ . Once C i has been found, equation (37) can be employed to 
find the initial value of (A j ) 0 - The remaining unknown initial condition, (A 2 ) 0 , follows imme- 
diately from equation (21) whether or not the previous two initial conditions have been found. 

The problem of determining the optimal trajectory has now been reduced to the isolation of 
the initial values of f3 y and f3 y and some considerations will now be directed toward their 
relationships to the Lagrange multipliers. 

It should be kept in mind that in this option, the use of the thrust orientation angle (ft) and 
its first two derivatives is something of an artifice, and is used only to derive values of the 
multipliers at the initial point. Indeed, once these values have been found it is no longer neces- 
sary to obtain values of fi or its derivatives unless we desire to do so for physical considera- 
tions.* 

To cement ideas, let us consider a concrete example wherein a circular orbit of prespecified 
altitude is desired. This may be assured if circular energy, zero radial velocity, and altitude 
are all matched at the end point. Thus, we have only three end conditions that must be met, 
while at the initial point we have four degrees of freedom — C i , (Aj )o , (A 2 )o , and burning time. 
These first three correspond to initial choices of f3 y /3, and Thus, there is one degree of 
ireedom available, and we might fix this by any of several possibilities. We may for instance, 
require that /3be continuous across staging. We could also determine the initial value of /3 such 
that the angle of attack at cutoff is zero or some other prespecified value. Finally, we could 
determine the initial value of jS such that our payoff function (payload, for instance) is a maxi- 
mum with respect to this value, etc. In this case the necessary end point constraints associated 
with the Mayer problem are applicable to the determination of j3 at the initial point, or we could 
apply a numerical isolation procedure. 

The situation is slightly more complicated when we consider the isolation of C t and (Ai) 0 
by approximating techniques with equations (34) and (37), respectively. Once f3 0 and (A i) 0 
have been chosen, all terms of equation (34) are known (from booster cutoff conditions) except 
j8 0 > j6 0 , and C 1 . /3 0 and are inserted from the approximate turn program and the value of C x 

determined. It should be noted, however, that this equation is a quadratic in ^S 0 and we can thus 
have two values of /3q that yield the same C x . Furthermore, the value of /3 0 could, at this point, 
be arbitrarily set and a proper choice of Cj still found by variation of j3o . On the other hand, 
(Ai ) 0 is a linear function of /3 and C x , but is independent of j3 . It would seem that we are not 
immediately guaranteed a unique pair [fi 0 , /3 0 ] for a given pair [ ( A 2 ) 0 , C j ] . 

It is interesting to note that if (A j ) 0 is arbitrary (i.e., only two end conditions are of in- 
terest), then there are two values of /3 which give equivalent end results once the variational 
trajectory is complete. These two values may not both be physically realistic. For example, 
several cases were found where the value of was negative, presumably indicating that a 
retrograde orbit** was ahead. Nevertheless, normal orbits were attained just as if had been 
initially positive. 

A few further comments are in order. The first of these is that if a multi-stage vehicle is 


* Once the variational turn program has been initiated, equations (12), (13), (21), (24), (29), 
and (38) are integrated numerically. Equation (21) may be used to eliminate /3 from equations 
(12), (13), and (29). From this time point, we may then use equation (21) to determine jS, equa- 
(37) to determine /3, and equation (34) to find at any point along the trajectory. The angle of 
attack, a, may be found by use of equations (5) and (11) if it is of interest. 

** Assuming that the vehicle was initially fired with the planetary rotation. 
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considered we have violated the initial assumptions used in the formulation of the Mayer pro- 
blem since there now exist discontinuities in the mass and thrust of the vehicle which introduce 
discontinuities into the Lagrange multipliers. The interested reader is referred to Reference 4. 

Secondly, insofar as magnitudes of the multipliers are concerned, a few general statements 
can be made. The signs of (Xi ) 0 an d C x are negative for orbits that are with the Earth's rota- 
tion while (X 2 ) 0 positive. The magnitude of (Ai )o is usually between - 0.1 and - 0.5, but 
the actual value is critical. The magnitude of C 1 is generally bounded by — 10 5 < C i< — 10 6 .* 
The values of /3 0 , j8 0 , and /3 0 are far more easily available from an intuitive viewpoint. For 
example, it is evident that a rate of 5°/sec (or even l°/sec) is extremely high, whereas a value 
of 50° is not unreasonable for /3 at booster cutoff. 

It may be seen that most realistic end conditions can be met without the assumption of 
coasting between stages. This is often not the case with other turn programs. The gravity turn, 
for instance, is shown to require such coast periods in Reference 3* It is also of interest to 
note that it is shown in this source that better performance results when no coast periods occur. 

Finally, it is usually possible to optimize the cutoff conditions of a non-optimal booster 
program with respect to this turn program. Such a procedure has been found to yield significant 
performance increases. Such an optimization is also possible in the alternative option treated 
below. 


OPTION TWO 


Another possible solution to equation (34) (which is actually a special solution of Option 
One^is found by setting C x = 0 and assuming that (0*0 for all values of t. For this case 
we have 


CiM • T ■* • • * 

— j- sin 2 - 4 r 0 2 + — sin/3 tan ft + r (jS tan /3-2/3 2 -6/3i/0 


■f 2 r fi tan /3 = 0 (39) 

* Equations (12), (13), and (39) are sufficient for the complete solution of the ascent problem. 
We may solve equation (39) for /3 and simultaneously integrate these three equations which do 
not involve any Lagrange multipliers . 

It should he remembered that equation (39) is a special solution of the optimal turn program, 
but it is not an approximate solution except insofar as the basic assumptions limit all solutions. 

The problem of discontinuities in the multipliers across staging does not arise in this 
option. 

These advantages are partially offset by the loss of one degree of freedom. We now have 
only /3 0 , $0 » and the burning time at our disposal since j§ 0 is specified by equation (39). 
These three values are sufficient to attain a circular orbit of specified altitude, but a free 
choice /S 0 niay not always be realistic. In this case we can perform certain other tricks to ob- 
tain more than two end conditions. We could shape the booster trajectory such that no discon- 
tinuity need occur in the thrust orientation angle, or use coast periods of unspecified lengths 
between one or more stages to accomplish the desired ends. 


* C x may range down to zero (see Option Two), but die probability of obtaining values be- 
tween - 10 5 and 0 is small. 
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Equation (39) may have a closed form solution, but none has been discovered thus far. The 
solution 

P -t* (40) 

where is a constant of integration, is irritatingly close to a general solution of this equation. 
It is a solution only for the limiting case when the mass of the primary body vanishes (M = 0); 
in this limit, however the physical meaning of the solution becomes questionable. 

A recent derivation of Leitmann (Ref. 5) obtained a solution to a quite different problem 
which is similar in form to equation (40). This solution, based originally on Reference 2, was 
first obtained by Lurie (Ref. 6). It is interesting to note that Leitmann’ s solution uses a tech- 
nique similar to that employed in the derivation of equation (34).* 

As a final comment, one often finds variational treatments of optimal trajectories which are 
formulated using only equations (1), (2), and (3). That is, the ,f range” constraint is deliberate- 
ly excluded. While the Lagrange multipliers of this report do not correspond (directly) with those 
used in such formulations. Option Two is equally as general and gives a much more useful re- 
sulting form. 


OTHER OPTIONS 


Equation (34) may have other possible solutions. The first such possibility that would seem 
to be of interest is obtained by setting X x (t) = 0, and requiring that C x ^ 0. Equation (34) may 
then be satisfied if equation (40) is fulfilled. Further probing into the basic set of equations of 
motion and the necessary Euler conditions indicates that further assumptions are necessary. The 
result of these further assumptions is that only unth rusted radial decent can satisfy all restric- 
tions. 

If Aj (t) = 0 and Cj = 0 we lose all information about the steering program. 

No other solutions to equation (34) have been found at the present time. 

SECTION VI. AN APPROXIMATE TURN PROGRAM 

The usefulness of Option One is a function of how accurately the optimal turn program can be 
approximated. This problem is, by its very nature^ rather iterative and the following material is 
simply a fair first guess. 

Ve shall assume that some trajectory analysis has given an approximation of the time to 
orbit and the maximum angle of attack that was necessary to obtain orbit. Using these data we 
now “curve fit” a symmetric parabola to die data as follows: 

Assume that a (*) is given by 

a = K i °~2 ^ (41) 

where K x and K 2 are constants, t 0 and t f are the initial and final times over which the turn 
program is applied, and t is the time along the trajectory. We now proceed to dete rmine K t and 

k 2 . 


The author was unaware of this reference during the development of the present theory. 
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Assuming that a = a max at t - 


1 o + 1 f 


, we 


find 


^2 fl mai 


and, now, 


a =a mBI + Ki ( * - 


t o + j f 


If the program is initiated with a = 0 at t = f 0 ,* as is often the case, we have 


a + K 

max 


.(PP -0 


Thus 


Finally 


-«max ( t( _ 1q 


This approximate turn immediately requires that , since a 0 iS 2ero > and 

/3o = cto + $o 

The first two derivatives of equation (46), at the initial point, are 

* i a ^max 

1 1= t 0 t { - t o 


« . o max 

a = “ o — — p 

f — *0 ( t { ~ *0 ) 


(42) 

(43) 

(44) 

(45) 

(46) 

(47) 

(48) 

(49) 


The initial values of fi 0 and /3 0 may now be evaluated from 


/3 0 - «o + *o 


and 


^0 =a o + *0 

The value of & 0 , from equations (2), (5), and (6), recalling that a 0 = 0, is 
* 


GM - r 0 fro + t-q <A q ) 

r iJVTiJi ) 


(50) 


(51) 


(52) 


* This also forces a to be zero at t = t { by symmetry. 
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By differentiation of equation (2) and similar substitutions, we find 


a 


/T5 

V r 0 + 


r 2 Jp. 

r o Vo 






CM 




r2 o r 0 ) 


« * f 

2 r Q xjj 0 

GM 

2 

GM 

+ 1 

'• l 

Jo ( T 'o+ T 0 t i'o> - 


u r o( r “o 2+ r 0 •A’o 2 ' 1 J 


( 53 ) 


This parabolic approximation could be supplanted by first employing Option Two and ob- 
taining a solution for Option One, with the results thus obtained. In this case, C i would be 
changed to a non-zero value with a simultaneous modification of a 0 or coast time toward their 
desired values. 

A gravity turn approximation (Ref. 3) may also be substituted for the parabolic function 
given above. 


SECTION VII. CONCLUSIONS 

The foregoing development has led to some rather interesting conclusions about optimal pla- 
nar trajectory programming of a constant thrust vehicle in vacuo. As pointed out in the intro- 
duction, these assumptions are usually employed in numerical studies of optimal upper stage 
trajectories. 

Option One is the most generally applicable case and the approximate turn program of Section 
VI allows a rapid isolation of the initial values of the Lagrange multipliers that usually bottle- 
neck optimal turn program studies. 

Option Two is of interest from at least two points of view. First of all, it is an optimal turn 
program that involves no Lagrange multipliers; to the best of the author's knowledge it has not 
been demonstrated previously. Secondly, it is probably simple enough to find application to 
adaptive guidance schemes for hardware vehicles. 

Both of the above solutions were investigated by numerical integration techniques. In all 
cases, they were found to agree with standard variationa 1 procedures. 
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